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Abstract. We consider the interaction of two vortex patches (elhptic 
Kirchhoff vortices) which move in an unbounded volume of an ideal incom- 
pressible fluid. A moment second-order model is used to describe the inter- 
action. The case of integrability of a Kirchhoff vortex and a point vortex 
by the variable separation method is qualitatively analyzed. A new case of 
integrability of two Kirchhoff vortices is found. A reduced form of equations 
for two Kirchhoff vortices is proposed and used to analyze their regular and 
chaotic behavior. 



1 Introduction. 



The simplest example of planar vortex motion of an ideal fluid, other 
than that described by the point vortex model, was suggested by Kirchhoff 
P] . He established that an elliptical vortex patch with semiaxes of a, b and a 
uniform vorticity u inside uniformly rotates around its center with an angular 
velocity of Q = {ujX)/{1 + A)^, A = a/b. In this rotation, the fluid particles 
are involved in the absolute motion with a double angular velocity (Lamb, 
1932). Lord Kelvin (1880) and Love (1893) showed that the Kirchhoff vortex 
is neutrally stable if and only if a/b < 3. 

S. A. Chaplygin in 1899 generalized Kirchhoff 's solution by introducing 
a uniform background vorticity into the unbounded fluid surrounding the 
elliptic vortex (this is the so-called Couette simple shearing motion). He 
established that the vortex will rotate with a certain angular velocity and 
change its contour (pulsate) in accordance with a certain law, which he ob- 
tained by integrating a system of two nonlinear differential equations. He 
also analyzed in detail the behavior of pressure in the fluid as a function of 
time. 

Kida [2| and Neu 0| generalized Chaplygin's solution while being unaware 
of his work |3]. Chaplygin's works in many fields of mechanics are poorly 
known outside the Russian speaking world. Forgotten Chaplygin's works on 
two-dimensional vortex structures are discussed in the recent historical review 
by Meleshko and van Heijst A superposition of the solutions obtained by 
Kida and Neu is presented in the book by Newton [0]. It was found that the 
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dynamics of a Kirchhoff vortex can be reduced to a one-degree-of-freedom 
Hamiltonian system if the velocity of the external flow can be described by 

U = Ui + U2 + U3, 

ui = {-fix, -72?/, 0), 

(1) 

U2 = (0,0,732;), 

U3 = (-741/, 74a:, 0) 

under the condition 71—72 + 73 = derived from the incompressibility condi- 
tion, and. It is assumed that the elliptic patch is the section of the elliptical 
cylinder determined by the plane xy. Equation ((H) yields Chaplygin's solu- 
tion when 7i = 72 = 74 7^ 0, 73 = 0; Kida's solution when 73 = 0, 71 = 72; 
and Neu's solution when 74 = 0. Note that Ui corresponds to the uniform 
deformation field (which is induced at a distance by one point vortex 
jH]); the field U2 corresponds to the extension along z axis, and the field U3 
corresponds to uniform background rotation. The Hamiltonian system that 
describes this case has the form 

(2) 

where the dot denotes differentiation with respect to the new time r = 
u;(t)(AV(A^ - while u{t) = uoe^^^-f^^^ and ujq is the initial vorticity 
in the Kirchhoff vortex. If 71 = 72, the system (j2|) is autonomous and can be 
explicitly integrated. Qualitative analysis of this system is given in [2], |S!, 

When the coefficients in are functions of time (for example, period- 
ical), we obtain (at 71 = 72) a Hamiltonian system with one and a half 



degrees of freedom with periodic perturbations, which was studied from the 
viewpoint of sphtting of separatrices and the appearance of stochastic behav- 
ior [Oj , JD] • The problem of advection of a passive particle of the fluid in the 
Kida flow was discussed in yjj. Examination of Poincare sections allows us 
to conclude that the advection is chaotic, i.e., Lagrangian turbulence exists. 
The stability of the elliptic vortex in a Kida flow (uniform deformation) was 
studied in Ref. [12]. 

Finally, we will mention the generalized Kirchhoff solution obtained by 
Polvani and Flierl ^Hj and corresponding to a family of embedded confocal 
elliptic vortices with the appropriate distribution of vorticity. In that work, 
the stability conditions, which generalize the Kelvin-Love conditions, were 
obtained for a system of two confocal elliptic vortices. 

2 Moment model of Kirchhoff vortex inter- 
action (vortex patch dynamics). 

The second-order moment model ^1] represents a higher level of approxi- 
mation in the description of hydrodynamic vorticity as compared with the 
point vortex model. This model considers elliptic patches with a specified 
vorticity, which move in a two-dimensional incompressible fluid. This the- 
ory was suggested by Melander, Zabusky, and Stychek (MZS model) in Ref. 
|14j . Rotation of two Kirchhoff vortices in the presence of central symmetry 
is analyzed in . The nonsymmetrical situation of interaction between two 
vortex patches is discussed in ^H], where the same moment model, which 
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describes merging of vortices, is compared with a pseudospectral model (in 
which Euler equations with a weak dissipation are solved). 

Let us describe this model in more detail. The second-order moment 
model, describing Kirchhoff vortex interactions, is derived from two basic 
assumptions: 

1. the distance between vortices in the process of evolution is much greater 
than the vortex size, therefore the vortices will maintain their elliptic 
shape; 

2. the third and higher moments in the expansion of the Hamiltonian can 
be neglected. 

Under these assumptions, the equations of motion of elliptic vortices can 
be written in the Hamiltonian form with nonlinear Poisson brackets 



Xk = {xk, H}, ijk = {vk, H}, (pk = {(pk, H}, Afc = {Afc, H}, 

(a^i, Vj} = p-dij, Wi, Xj} = Y-S^i-X^ ^'^' 
and the Hamiltonian is 



(3) 



H = Hi + H2 + -f^s, 

H. = -ltni-^^^ i/. = -^f'r.r,inM.„ 



(4) 



k=l k,p 

1 — A^ 

+ cos(2(^fcp - ^k)) 

■^k 

where F^, Sk are the total intensity and the area of the elliptic vortex with 
index k; Mkp is the squared distance between the centers of the k-th and p-th 



vortices {M^p = {xk — XpY + {y^ — Up)^)] is the slope of the fc-th elhpse 
with respect to the and 6kp is the angle between the x-axis and the 

straight line connecting the centers of the fc-th and p-th ellipses (see Fig. 1). 

The area of each ellipse remains constant {S^ = const) by virtue of 
Kelvin's theorem on the conservation of circulation in an ideal medium [T7] . 
jl8j : therefore, these areas are parameters of the model considered. 

The components of the Hamiltonian have the following physical meaning: 

Hi describes the action of an elliptic vortex on itself, see 0; 

H2 describes the interaction between equivalent point vortices; 

describes the interaction between vortices with second-order moments 
taken into account. 

In addition to the Hamiltonian H, equations (jHI) have a noncommutative 
set of first integrals 

N N N 

Q = ^ TfcXfc, P = ^ '^kVk, / = ^ Tfc 

k k k 

which correspond to the translational and rotational invariance of the system 
in the absolute space. 

Integrals Q, P, I satisfy the following commutation relations: 

{Q,P} = 5^7., {P,I} = -2Q, {Q,I} = 2P. (6) 

i=l 

Therefore these relations are not sufficient even for the integrability of a 
system of two Kirchhoff vortices; this system will be further reduced to a 
system with two degrees of freedom. However, the problem of the dynamics 
of a Kirchhoff vortex and a point vortex (a system with three degrees of 
freedom) is integrable. 
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2 , 2 , Ski + 



Xk + yk + 



47r Xi 



, (5) 



A more general model of interaction between elliptic vortex patches was 
suggested by D. Dritscliel and B. Legras ^Hj, [20] • In some sense, this is an 
intermediate model between the MZS-model and an exact description deter- 
mined by the contour dynamics method (in the contour dynamics method, 
the elliptic shape of the patch does not persist, while the model [ini, [201 
is derived with no allowance made for the velocity field responsible for the 
nonelliptic part of the interaction). However, this model is more complex 
and in many cases it is sufficient to use the MZS-model. 



3 Interaction of a KirchhofF vortex with N 
point vortices. Integrable case for = 1. 

Let us denote the coordinates of the center of the elliptic vortex by (xo,i/o) 
and those of the point vortices, by {xi,yi), i = 1,. . . ,N] the equations de- 
scribing the dynamics of this system also can be written in Hamiltonian form 
with a Poisson bracket and the Hamiltonian in the form 



{xi.yi} = T-^5ij, i,j = 0, 1, . . . , A^, {(p,X} 



in 



A2 



TqS 1 - A2 
H = Hi + H2 + H3, 



1 n I 1 ^ / 

k,p=0 

H-^ = ^ ; COS 2i6k — 

^ 167r2 Mfeo A ^ 

k=l 

where Fq is the intensity of the Kirchhoff vortex with the semiaxes ratio 
of A with the angle 9?, determining its orientation (see Fig. 1), 5* is the 
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area of the ellipse, 6^ is the angle between the x-axis and the straight line 
connecting the center of the Kirchhoff vortex and the k-th point vortex, and 
Mkp = (xfc — Xp)^ + {yk — ypY- Hereafter, we will assume without loss of 
generality that the intensity of the Kirchhoff vortex is positive (Fo > 0). 

The integrals of motion corresponding to the group E{2) of motions of 
the plane are determined by relationships (j^I), where it should be assumed 
that 5*0 = S*, and = 0, i = 1, . . . , A^; their commutation relationships are 
analogous to ©• As a consequence of the existence of integrals, we obtain 

m 

Proposition 1 The system of an interacting Kirchhoff vortex and one point 
vortex {N = 1) is completely integrable. 

This was first shown by Lebedev and somewhat later was established 
independently by Riccardi and Piva [22]. Here we will give a geometric 
analysis of the motion of vortices to improve the results of fH] . 

For explicit integration and qualitative analysis, we perform a reduction 
to one degree of freedom. Let us consider new (relative) variables 

^ = 2(^-¥p), p = 1c(a + ^), z = M,o = ixi-xof + {yi-yof, (8) 

where c = Fo5'/87r. All these functions commute with the integrals (0), 
i.e., they are invariants of the group of motions of the plane E[2) and are 
closed with respect to the Poisson bracket 

{^,p} = i, {^,z} = -A{r,' + r^'), {p,z} = o. (9) 

The Poisson structure Q has a linear Kazimir function 

D = Toz + A{l + a)p, a = Fo'Fi. (10) 
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Eliminating z with the use of equation (fTUI) . we obtain the Hamiltonian 
of the system (up to the constant) with one degree of freedom: 



p2 / I 2 _ ^2 

H^ = ^\ -ln(c + p) - 2a ln(D - 4(1 + a)p) ^ 80—^^7- cos ^ 

svr \ — 4(1 + (x)p 

(11) 

where, according to Q, the variables ■?/', p are canonical. 

The different signs in the Hamiltonian (jllll appear due to the nonunique- 
ness of the inverse transformation (jH)) for p(A), which has the form 



0<A<1 

A=<( 'r^2 (12) 

P+vVHZ, A>1 
c 

(the upper sign in (fTTj) corresponds to the case A < 1). 

The reduction by (to two degrees of freedom) for arbitrary can be made 
in a similar manner. The new variables for the reduced system can be chosen 
as follows: 

T/^fc = 2(4-(^), p=lc(A + A-^), 
Mn, = [xi-Xkf ^{Vi-ykf, i = 0,...,Ar, k = \,...,N. 

Qualitative analysis of the relative motion for = 1. According to 
((HI), the domains of the variables ^z^, p are determined by the inequalities 

< < 47r, c < p. 

Each point of this half-strip corresponds to a pair of possible mutual 
arrangements of the elliptic and point vortices (fT^ . corresponding to A < 1 
and 1 < A (see Fig. 2). 
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It can be shown that the relation 



p = c, i. e. X = 1, 

corresponds to the case when elliptic vortex becomes circular (Rankin vortex ). 

The trajectories of the reduced system are determined by the contour 
hues of the Hamiltonian (lllj): therefore, the trajectories of the system in 
domains A < 1 and A > 1 can be obtained with the help of the substitution 

^p ^ 'ijj + n, p ^ p- 

Moreover, the Hamiltonian (jllll is 27r-periodic with respect to ip. Therefore, 
we will restrict ourselves to the analysis of the contour lines of the Hamilto- 
nian H_ in the domain 

< ^ < 27r, c < p. 

Normalizing the variable and the integral 

p = cy, D = Ac\l + a\D (13) 

and eliminating "insignificant" constants in ()11|). we find that the trajectories 
of the system are determined by the level curves of the function 

~ ~ 2(y — 1 

E_ = - ln(l + - 2a ln(D -y) ^ cos^, 1 + a > 0, 

1 + a D-y 

~ ~ 2a - 1 

E_ = - ln(l + w) - 2a ln(D + w) - cos V^, 1 + a < 0. 

|l + a| B^y 

(14) 

According to [14], the equations describing the dynamics of a Kirchhoff 
vortex interacting with a point vortex are valid only at a sufficiently large 
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distance from the Kirchhoff vortex. Without specifying the exact domain 
where the results thus obtained are apphcable, we give here a trajectory of 
the reduced system flllj) and mark the domain occupied by the Kirchhoff 
vortex. In accordance with (fTUI) . (fT^ . the elhptic vortex in the plane ip, y 
occupies the domain determined by the relationships 
^ 2 

D-y<— — {y- ^J%p - IcosV'), 1 + a > 0, 

D + y<- -{y - Jy'^ - 1 cost/)), 1 + a < 0. 

|1 + a| 

The domain occupied by the Kirchhoff vortex is shaded in the figures below. 



Stability of the circular vortex. Let us make the canonical change of vari- 
ables 

y = 1 -\ , ip = arctan — . (16) 

Now the Hamiltonian H14|) near y = 1 can be written as 
~ 2V2a ib-l-Aa, o 

i7 = const ^ V ^ iu^ + v^) + ..., 1 + a > 0: 

{l + a){D-l) 2 D-1 

~ 2V2a ib + l + Aa, o 

if = const H V ^ iu^ + v^) + ..., l + a<0. 

{l + a){D + l) 2 + i ' 



Thus, after the identification of (|16|). we obtain that the origin of coordinates 
u = V = \s not a fixed point, i.e., in the presence of point vortex, the circular 
vortex is locally unstable with respect to elliptic deformations. 

Let us qualitatively describe the structure of the phase portrait on the 
if), y plane for different values of the parameters a, D. As follows from (fTIHl 
and (HH), the phase space of the reduced system and the respective phase 
portraits change depending on the sign of (1 + a). Let us consider each case 
separately. 
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1 + a > (Fig. 3). In this case, tlie ptiase space is the rectangle on 
the plane ip^y: 

0<^<27r, l<y<b. (17) 

The opposite sides of this rectangle if) = Q and = 2tt are identified with one 
another. The segment y = 1 corresponds to the case of a circular Kirchhoff 
vortex, and y = D is the case when the point vortex is in the center of the 
elliptic vortex. In both cases the definition of the angle if) makes no sense 
(i.e., this is a singularity similar to the origin for polar coordinates). In this 
case, the phase space can be identified with a two-dimensional sphere. 

Remark 1 The identification for the sphere + ■Cl + 'Ci = 1 embedded 
into can be explicitly specified, for example, in the following manner 

6= ^1 = v/l-e|cosV^, 6 = ^1-^1 sin V^. (18) 

Fig. 3 gives a characteristic appearance of phase portraits corresponding 
to the three different parameter domains. The appearance of the portrait is 
completely determined by the critical points of the Hamilton function (fT^ . 
which, as follows from (fT^ . lie on the straight lines = 0, ip = n. In this 
case the point vortex lies on the principal axes of the elliptic Kirchhoff vortex. 

Each straight line ip = and = tt {in the nondegenerate case) contains 
either or 2 critical points of the function ()14|1 . 

The domain (fTTjl can contain either 2 or (nondegenerate) critical points 
of the function ()14|) that do not lie on the segments y = 1, y = D. Both 
critical points lie either on the straight line x = or on the straight line 
X = TT. The critical point nearest to the segment y = 1 always corresponds 
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to a stable fixed point of the reduced system. 

A criterion for the global stability of the Kirchhoff vortex at 1 + a > 
can be formulated within the framework of the model considered: 

if the Hamilton function has critical points in the domain (fTTj) . the per- 
turbations of the Rankin vortex remain constrained. 

Indeed, in this case, as can be seen from Fig. 3, there always exists an 
invariant curve that limits the extent of the Kirchhoff vortex deformation 
and inhibits its merging with the point vortex. 

1 + a < (Fig. 4). In this case, the domain of motion in the ^-plane 
is noncompact: 

<■?/'< 27r, y> max ^1, = y^. 

Carrying out the (canonical) change of coordinates 

y = yo^ 7^ — , V = arctan -, 

2 V 

we find that the phase space is the (m, f )-plane and the straight line y = yo 
corresponds to the origin of coordinates. The structure of the phase portrait 
(the shape of the domain occupied by the Kirchhoff vortex) depends on the 
sign of D. 

When D > 0, there is a stable periodic solution on the straight line if) = tt 
(Fig. 4a). For D < 0, an unstable periodic solution exists aX ip = Q (Fig. 46). 

l+a = (the case of a vortex pair). The transformation (fTT^ cannot 
be made in this case, therefore we carry out the normalization 

p = cy, D = Acb = FoMio = const. (19) 
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Thus, in this case, the distance between vortices remains constant (Mio = 
const), and changes take place only in the mutual arrangement of the vortices. 
The variables ip, y are determined only in the half-strip 

0<tp <27T, 1 <y. 

In this case the trajectories are determined by the contour lines of the Hamil- 
tonian if_ in (|TT|l . which, after the elimination of constant values can be 
written as 

= -ln(l + y) + iVP- Icost/-, (20) 

The domain occupied by the elliptic vortex on the domain of the variables is 
determined by the inequality 

D>2{y- VP _ icos^). (21) 

Whence it follows that for D < the entire plane (V',y) is occupied by the 
elliptic vortex; therefore, we assume that D > 0. 

For D < = \/22 + 10\/5, the phase portrait contains no periodic 
solutions (Fig. 5a). 

When D > D^, two periodic solutions appear in the phase portrait on 
the axis ijj = 0; one of these solutions is stable, while the other is unstable 
(Fig. 5b, c) 

4 Interaction between two KirchhofF vortices 

Reduction to a system with two degrees of freedom. The dy- 
namics of two KirchhofF vortices can be described by the Hamiltonian 
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which can be represented in the form 



H = Hi + H2 + H3, 



Hi 




H2 = -^^InM, 



riT2 



Sijl-Xj 
Ai 




cos 2 



{9 - ^2) 



M={xi-X2f+iy,-y2f, 6 = 9 



7^+021 = arctan 



1/2 -yi 



X2-X1 



(22) 



The particular solution of the system (j22I), at which 71 = 72 and the 
vortices are centrally symmetrical with respect to one another is suggested 
and studied by Melander et al. ^3] . In that paper, the conditions for merging 
of two vortex patches are given. These conditions, within the framework 
of the moment model, were found to be equivalent to the collapse of two 
Kirchhoff vortices, during which their centers coincide after a finite time. 

Let us consider the system of relative variables 



These variables commute with the integrals (jni); they are closed with 
respect to the Poisson brackets (jH)), and their commutation relationships 
have the form 



{V^,,p,} = 5,„ {V'i,4 = -(rr^ + r2-i), {a,z} = 0, z,j = 1,2. (24) 



The Poisson structure pH) has a linear Kazimir function (integral of 




(23) 



z = -M. 
4 
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motion) 

D = ^+(rr^ + r2^)(pi + P2). (25) 

As follows from (j21|), (f^ - in the case of the vortex pair (Fi + r2 = 0), 
the distances between the vortex centers remain constant. 

Eliminating z using the integral and expressing the Hamiltonian (j^^ 
in terms of the variables pi, we obtain the Hamiltonian of the reduced 
canonical system with two degrees of freedom 

if = --i ln(ci + pi) - / ln(c2 + P2) 

^^^Mn (D-(rr^ + r2-^)(pi + p2)) 



47r 



rir2 i^^^^. , ^MEA , 

' cos ipi =F -^^ — cos %)2 



4^p_(r-i + r2-^)(p, + P2)) \^ Ti 

(26) 

where the upper sign corresponds to the condition Aj > 1, and the lower sign 
to the condition < Aj < 1 (see. (fT^ ). 

Absolute motion. Let us find the equations determining the positions 
(xi,?/i) and {x2,y2) and orientation (pi, of two Kirchhoff vortices in sta- 
tionary space if their relative arrangement determined by the system fj26|) 
is assumed to be known: pi = pi{t), ipi = ip{t), i = 1,2. Direct calcula- 
tions show that the slopes of the principal (major) semiaxes of the ellipse are 
determined by the quadratures 
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If Fi + r2 7^ 0, the positions of the vortex centers can be found from the 
system of hnear equations with coefficients exphcitly depending on time: 



= — ^2]vf2 I (AxsmV'i + AycosV'i) 



^ / 2 ^ 



4co 



" n'^M'^ \ 4^^ (Ax cos ^1 - Aysmipi^ 



- ( Ax cos '?/^2 — A?/ sm '1/^2 



4C2 



:28i 



where Ax = Xi — X2, Ay = yi — y2, and Q, P are integrals (0). 



Integrable case Fi + r2 = 0. In the case of the vortex pair Fi = — F2, the 
Hamiltonian can be split into two independent Hamiltonians 

and the system can be integrated by the method of separation of variables. 
Thus, we obtain a new nontrivial integrable case for the vortex dynamics. 
This case of integrability was first mentioned in our review |2S!- As follows 
from equations (|^. in this case 



Xi — X2 = const, yi — y2 = const, 

and hence 

9 = arctan — — — = const. 

X2 - Xi 
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Thus, the length and orientation of the segment connecting the centers of 
vortices remain unchanged during the motion (see Fig. 1). 

Poincare sections and integrability. The system p6|l is not integrable 
in the general case. This is demonstrated by the chaotic trajectories con- 
structed at Fi = F2 with the use of a Poincare section shown in Fig. 6. It 
should be mentioned that there are few returning trajectories in the case of 
interaction between Kirchhoff vortices, which hampers the numerical analy- 
sis. Fig. 7 gives numerically plotted separatrices of a hyperbolic fixed point 
of the Poincare mapping. Their transversal intersection (established by com- 
puter analysis) is an indication that the problem of Kirchhoff vortex motion 
is nonintegrable. This fact was mentioned as a supposition in No an- 
alytical proof of the nonintegrability of two Kirchhoff vortices within the 
framework of the second-order moment model has been obtained yet. 
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Figure 1: Elliptic vortex patches 





A > 1 



A< 1 



Figure 2: Mutual arrangement of the elliptic and point vortices with respect 
to a chosen point on the contour of the elliptic vortex in the case where A < 1 
and A > 1. 
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e) a = 10, L>=1.5 



Figure 3: Phase portraits for l+a > (the domain occupied by the Kirchhoff 
vortex is colored gray). Figure 3d is an enlarged part of figure 3c near the 
lower segment. 
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Figure 4: Phase portraits for 1+a < (the domain occupied by the Kirchhoff 
vortex is colored gray). 
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Figure 5: Phase portrait for 1 + a = (the domain occupied by the Kirchhoff 
vortex is colored gray). Figure 5c is an enlarged part of figure 56 near the 
lower segment. 
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H= -0.238 



H= -0.24 



Figure 6: Poincare mapping for the reduced system (|^?)|) in the problem of 
two elhptic vortices. The section by the plane ipi = was chosen for the 
following parameter values: 7i = 72 = 1, Si = S2 = 0.3 and D = 22. 
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